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ABSTRACT 

We examine metal and entropy content in galaxy groups having Tx ~ 0.5 — 2 keV in 
cosmological hydrodynamic simulations. Our simulations include a well-constrained 
prescription for galactic outflows following momentum-driven wind scalings, and a 
sophisticated chemical evolution model. Our simulation with no outflows reproduces 
observed iron abundances in X-ray emitting gas, but the oxygen abundance is too low; 
including outflows yields iron and oxygen abundances in good agreement with data. 
X-ray measures of [0/Fe] primarily reflect metal distribution mechanisms into hot 
gas, not the ratio of Type la to Type II supernovae within the group. Iron abundance 
increases by x2 from z ~ 1 ^ independent of group size, consistent with that 
seen in clusters, while [0/Fe] drops by ~ 30%. Core entropy versus temperature is 
elevated over self-similar predictions regardless of outflows due to radiative cooling 
removing low-entropy gas, but outflows provide an additional entropy contribution 
below 1 keV. This results in a noticeable break in the Lx — Tx relation below ^ 1 keV, 
as observed. Entropy at -R500 is also in good agreement with data, and is unaffected 
by outflows. Importantly, outflows serve to reduce the stellar content of groups to 
observed levels. Specific energy injection from outflows drops with group mass, and 
exceeds the thermal energy for < 0.5 keV systems. Radial profiles from simulations 
are in broad agreement with observations, but there remain non-trivial discrepancies 
that may reflect an excess of late-time star formation in central group galaxies in our 
simulations. Our model with outflows suggests a connection between physical processes 
of galaxy formation and both pre-heating and enrichment in intragroup gas, though 
more definitive conclusions must await a model that simultaneously suppresses cooling 
flows as observed. 

Key words: galaxies: formation, X-rays: galaxies: clusters, galaxies: clusters: general, 
galaxies: abundances, methods: N-body simulations 



1 INTRODUCTION 

Galaxy clusters are the largest bound structures in the Uni- 
verse. This oft-repeated phrase captures both their extreme- 
ness as well as their simplicity. Clusters are virialized and 
nearly spherical, their baryons are mostly in the form of 
a hot intracluster medium (ICM), and their galaxies are 
nearly all early-types with old stellar populations. The ICM 
is close to hydrostatic equilibrium with a temperature near 
the halo's virial temperature, and is enriched to around one- 
third solar independent of cluster mass. The only significant 
subdivision of clusters comes from examining their central 
gas properties, where they split into cool core (or cooling 
flow) and non-cool core clusters. Their smaller counterparts, 
galaxy groups (usually referring to systems with X-ray tem- 
peratures below 2 keV or so), share many of these character- 
istics. The simplest models regard groups as mini-clusters. 



i.e. bound spheres of hot gas in hydrostatic equilibrium, 
whose prope r ties should be self-similar with clusters (e.g. 
iKaiser Ill986l : iNavarro. Frenk. fc White Ill995l ). 

But a closer inspection of group and cluster properties 
reveals many puzzles and complexities. For instance, gas in 
hydrostatic equilibrium near the centers of groups should be 
sufficiently dense to ha ve coolin g times significantly shorter 
than a Hubble time ( iFabian., Nu lsen, & Canizarcs 198^), 
and therefore gas should be cooling rapidly onto the cen- 
tral galaxy. Yet X-ray observations show a distinct lack of 
gas at temperatures below roughly one-third of the system 
temperature, a nd little cold gas or star fo rmation is seen in 
the center (see [Peterson fc Fabian 1 l2006l . for review). Evi- 
dently, some source of energy or entropy is acting to pre- 
vent gas cooling in the core. Recently, observations of X-ray 
bubbles in the ICM apparently blown by AGN-driven jets 
from the central galaxy offer a possible way to balance cool- 
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ing llMcNamara et al I l2007l : iDunn fc Fabian 1 12008|) . Cos- 
mic r avs iIGuo fc Oh 1 12008|; ISiiacki et al. 1 12008|). conduc 



tion (|Kim fc Naravanl l2003l : |jubelga .s. Spri ngel. fc Dolag I 
20041'). and gravitational heating jKim fc Naravan I 20031 : 



Burns et al. 



20081 : iDekel fc Birnboim I l2008h may also be 



viable mechanisms. 

Another indication of non-gravitational energy injection 
into intracluster gas comes from scaling relations between 
cluster X-ray luminosity (Lx), X-ray luminosity- weighted 
temperature (Tx), and galaxy velocity dispersion (a). The 
self-similar model including only gravitational heating pre- 
dicts Lx oc Tx, which is confirmed by hydrodynamic sim- 
ulations that do not include ga s cooling or processes as- 
socia t ed with galaxy formation (iNavarro. Frenk. fc White I 
1 19951 : lEke. Navarro, fc Frenk1ll998ll. Observations ind icate 



Lx oc Tx for clusters l|White. Jones, fc Forman 1997 ), an d 
an even steeper relation in groups (e.g. Ixue fc Wu Il2000l ). 
This has led to speculation that some feedback process has 
added a significant amount of non-gravitational energy into 
the ICM, which would puff out the gas distribution and lower 
the X-ray lumi nosity. Such e nergy injection was dubbed 
"pre-heating" bv lKaiser 1 (|l99lh . Understanding the amount, 
epoch, and origin of pre-heating has been the subject of 
many investigations. 

An influential paper by iPonman. Cannon, fc Navarro I 

l|l999ll claimed that systems above 1-2 keV showed rising 
entropy with temperature, while below that temperature all 
groups had a similar entropy of ~ 100 keV cm^ . Here entropy 
is defined as the X-ray weighted temperature divided by the 
electron density at one-tenth the virial radius to the two- 
thirds power. While subsequent observations showed that 
this claimed entropy fioor is instead merely a more-gra dual- 
than-expected decline (e.g. iHelsdon fc Ponman I I2OOOI ). the 
question of what elevates the entropy of gas in smaller sys- 
tems remains unanswered. Entropy has now become the de 
facto quantity discussed when quantifying pre-heating. 

A third, more oblique indication of non-gravitational 
energy input into the ICM is the presence of metals. 
The amount of metals, roughly one-third solar in clus- 
ters regardless of size, is large when one considers that in 
the most massive clusters only a few percent of baryons 
have formed into stars. It is further remarkable that 
this enrichment level is similar in clusters of all masses 
despite significant trends with mass in the amount of 
stars formed (e.g. Gonzalez. Zaritskv. fc Zabludoff] l2007l : 
iBalogh et ai~l 120071 '). Enriching the ICM to these lev- 
els represents a non-trivial challenge, because stars 
in cluster galaxies can themselves o nly provide a 
metallicity of about one-sixth solar (jPortinari et al. I 
I2OO4I '). Accountin g for intracluster stars can a ll eviate 
this disrepancy (Zaritskv. Gonzalez, fc Zabludoff I |2004| : 



ISivanandam et al. I 2008[) , but the metals still must be re- 
moved from galaxies yet be retained in the hot ICM gas. 
One possibility is that gas strip ped off of infalling galaxie s 
can sufficiently enrich the ICM (|van Kampen et al. II2007I ): 
if so, the enrichment process would not be energetically im- 
portant. However, it is clear that the diffuse intergalactic 
medium (IGM ) is enriched at as early as z ~ 6 by pow- 
erful outfiows (lAguirre et al. II2OOII : lOppenh eimer fc Dave I 
I2OO6I '). so perhaps such outflows could also be responsible 
for enriching intracluster and intragroup gas. If so, outflows 
might simultaneously provide some level of pre-heating. 



In short there are three current puzzles in understand- 
ing the physics of the ICM: 

• Cooling flow problem: Cooling rates of hundreds to 
thousands of solar masses per year are expected based on 
simple cluster models, yet the measured amount of cooling 
gas is at least an order of magnitude less. 

• Entropy problem: Entropy levels are progressively more 
elevated over self-similar predictions towards smaller sys- 
tems. 

• Metallicity problem: It is unclear how metals are re- 
moved from galaxies in such a way as to enrich intracluster 
and intragroup gas as observed. 

Occam's razor leads one to attempt to solve these 
issues with a single physical mechanism. Unfortunately, 
this is difficult. For instance, the solution to the cooling 
ffow puzzle cannot invoke significant star formation (and 
hence metal production), because little recent star forma- 
tion is seen in cluster galaxies. Models of entropy injection 
from supernovae seem to suggest that it is insufficient to 
solve the entropy p r oblem llBalogh. Babul, fc Pattonl Il999l : 
iBower et all I2001I : iPipino et al 1 12002|) . These davs. it is 
fashionable to mostly neglect the metallicity problem, as- 
cribing it to a change in stellar yields or perhaps even 
the initial mass function (IMF), and attempt to solve the 
cooling ffow and entropy problems with a single energy 
source: heating from active galactic nuclei (AGN). However, 
this has difficul ties as well, because the entropy require - 
ments are large l|Babul et al. II2002I : [McCarthy et al. \\200i ). 
though more recently analytic models th at use AGN in a 
self-r egulated way have had some success (jMcCarthv et alH 
120081 ). 

In this paper we focus on the entropy and metallic- 
ity problems, leaving aside the cooling fiow issue for now. 
Previous theoretical studies of the entropy problem have 
generally assumed some amount of energy or entropy injec- 
tion at an early epoch, and then attempted to constrain the 
amount and epoch of pre-heating by evolving systems (either 
analytically or numerically) to compare w i th present-day 
clust ers (e.g. iBabul et al~fl2002l : iKav II2004I : iBorgani et ah] 
l2005h . Few attempts have been made to understand or 
model in detail possible sources of pre-heating. Observation- 
ally, there have been attempts to connect pre-heating and 
metal enrichment, but with limited success. iRenzini I (|200ll ) 
argued that galactic winds must be responsible for ICM 
enrichmen t, but they could only add a neglig i ble am ount 
of energy. IPonman. Sanderson, fc Finoguenov I l|2003l . here- 
after PSF03) found that neither isentropic pre-heating nor 
cooling models can reproduce the entropy structure of in- 
tragroup gas. Hence the physics of pre-heating remains a 
mystery, and the constraints on its amount and epoch re- 
main only loosely connected to physical process of galaxy 
formation. 

Here we take a different approach. We begin with a 
model for galactic outffows within hierarchical galaxy for- 
mation simulations that is constrained by non-ICM related 
observations, such as high-redshift IGM enrichment and the 
galaxy mass-metallicity relation. We then study how such 
outflows impact the enrichment and entropy level of intra- 
group gas at z = 0. We are limited to investigating galaxy 
groups and not clusters owing to our numerical capabili- 
ties, but this is actually desirable because groups fall in a 
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particularly interesting regime for pre-heating and enrich- 
ment processes. The obser ved turn-down at ~ 1 keV in the 
Lx — Tx scaling relations (iXue fc Wu Il2000l ) indicates that 
pre-heating energy input is comparable to the system's grav- 
itational energy at group scales. Groups have a signficantly 
larger stellar fraction, yet their metallicity appears to be 
comparable to or even lower than clusters, indicating that 
metal transport processes are crucial. Groups are more gen- 
erally important because they contain the majority of galax- 
ies in the Universe, and simple dynamical arguments suggest 
that group environments are the most effective for morpho- 
logically transforming galaxies through mergers. Hence there 
appears to be a lot of action at the group scale. Observations 
of intragroup gas are rapidly improving, and this provides 
important constraints on models of pre-heating and enrich- 
ment. Investigating such constraints is a central theme of 
this paper. 

We have recently identified a model for galactic out- 
flows that is remarkably successful at reproducing a wide 
range of observations of cos mic chemical enrichment. In 
lOppenheimer fc Dave ] l|2006l . hereafter OD06) we used cos- 
mological simulations that incorporate enriched kinetic feed- 
back to show that a certain class of outflow models are par- 
ticularly successful at matching detailed properties of IGM 
enrichment, as tra ced by C iv absorbers in quasar spectra 
from 2 ~ 6 ^ 2. In lFinlator fc Dave I (fioO?) we showed that 
the exact same model was concurrently and fairly uniquely 
successful at reprod ucing the observed mass-metallicity re - 
lation in galaxies. In lDave. Finlator fc Oppenheimerl (|2006l ) 
we showed that such winds are able to reproduce the lumi- 
nosity functions of early galaxies, significantly suppressing 
earl y star formation as need ed to solve the overcooling prob- 
lem. iBouwens et al. I |2003) showed that such a model best 
reproduces the early cosmic luminosity density evolution. 
The success of a single outflow model in matching such a 
wide range of observations, albeit mostly at z > 2, is com- 
pelling, and merits investigation into its impact on intra- 
group and intracluster gas. 

The outflow model we have found to be most success- 
ful follow s scalings predicted by a momentum- driven wind 
scenario (|Murrav. Quatert. fc Thompson 1 12005| ). In such a 
model, the wind speed scales with the galaxy's circular ve- 
locity, while the mass loading factor, i.e. the amount of mass 
ejected in an outflow as compared to the star formation rate, 
scales inversely with it. This yields a constant momentum in- 
put per unit star formation. We have found that the generic 
property making this outflow model successful is that early 
galaxies drive out an amount of mass comparable to or ex- 
ceeding that formed into stars, at speeds well above the 
escape velocity of their host halos. While this might seem 

suprising, it appears to be required for early IGM enrich- 

I — 1 ' \ 

ment JA guirrc ct al. 2001), global star formation suppres- 
sion (ISprin gcl & Hcrnquist 2003lJ), an d establishi ng the gas 
and metal content of high-z galaxies (|Erb I fioosi ) . Intrigu- 
ingly, this scaling of wind velocity with gal axy size is also 
observed in local starburst-driven outflows (jMartin I l2005l : 
iRupke. Veilleux fc Sanders ' '20051. and at 2 ~ 1 in star- 
forming galaxies IWeincr et al. 2008), providing a connec- 
tion between rare p resent-day outflows and the ubiquitous 
ones at high-z (e.g. ISteidel et al~ll2004l ). One notable issue 
is that for sufficient mass to be ejected, the energetic re- 
quirements are large, comparable to the available supernova 



energy l|Oppenheimer fc Dave |[200^ . Given expected radia- 
tive losses, it may be that another source of energy besides 
supernovae is required, such as photons from young stars in 
low-mass galaxies, or AGN in high-mass galaxies. The suc- 
cesses of these momentum-driven wind scalings suggest that 
we now have a plausible model, albeit phenomenological, of 
how outflows move mass, metals, and energy on large scales 
within a hierarchical structure formation scenario. 

In this paper we show that cosmological hydrodynamic 
simulations incorporating momentum-driven outflow scal- 
ings are able to simultaneously enrich intragroup media, 
reduce the stellar fraction to observed levels, and match 
the entropy levels of observed groups. We compare simula- 
tions with and without winds to assess the impact of outflow 
metal and energy injection. We find that radiative cooling is 
mainly responsible for breaking the self-similarity between 
groups and clusters, but that outflows provide an additional 
entropy boost that becomes important in the poorest groups 
{Tx ^ 1 keV). A key aspect of outflows is that they sup- 
press star formation, because cooling-only models can match 
intragroup medium entropy levels only at the expense of ex- 
cessive stellar content. Our results suggest that pre-heating 
and enrichment in groups today arise from the same out- 
flows that are required to distribute metals in the cosmos 
and regulate early star formation. 

In 321 we describe our Gadget-2 simulations with our 
chemical evolution and outflow modeling. In ^ we study 
the enrichment properties of intragroup gas in these models, 
including alpha element ratios and evolution from z = 1 ^ 
0. In S|4]we examine entropy levels and scaling relations in 
intragroup gas, and quantify energy injection from outflows. 
In 3S1 we compare radial profiles of physical properties to 
data, and study the impact of outflows on such profiles. We 
summarize in SJS] 



2 SIMULATIONS 



2.1 Simulations with outflows 



We run cosmolog ical hydr o dynam ic 
ing Gadget-2 l|Springel I |2005| ). 



simulations us- 
a cosmological 
Tree-Particle Mesh-Smoothed Particle Hydrody- 
namics cod e including radiative cooling and star 
formation (ISpringel fc Hernquist . .2003a. ) . Our ver- 



sion (|OppenIieimer fc Dave 20081. hereafter OD08) als o 



includes metal-line cooling l|Sutherland fc Dopita I Il993l ). 
uses a (s patially-uniform) c o smic photoionizing background 
given bv iHaardt fc Madau I (|200ll ). 

We include the effects of galactic outfiows and the cre- 
ation of heavy elements as describe d in OD08. This builds 
on our implementation detailed in Oppcnhei mer fc Dave I 
(2006), following Sprin gel fc Hernq uist (2003a), with three 
major improvements: (1) We track enrichment contributions 
from Type II supernovae. Type la supernovae, and mass and 
metal loss from AGB stars; (2) We individual follow car- 
bon, oxygen, silicon, and iron enrichment; (3) We use galaxy 
baryonic masses (rather than the local gravitational poten- 
tial) to estimate the galaxy parameters from which outflow 
properties are derived, using an on-the-fly galaxy finder. The 
upshot of these improvements is a more accurate tracking of 
metal evolution in various phases down to z = G, as well as a 
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more realistic implementation of the dependence of outflow 
properties with galaxy size. Our chemodynamical model is 
broadly similar to tha t used to study I CM enrichment using 
Gadget-2 bv iTornatore et al. I (|2007|) , thou gh some details 
differ. We employ iBruzual fc Charlotl (|2003l ) models with a 
Chabrier IMF for computing stellar evolution. 

The differences in our new implementation versus the 
one in OD06 are motivated primarily by a desire to treat the 
momentum-driven wind model more accurately. In our old 
model, we used the local gravitational potential as a proxy 
for galaxy velocity dispersion. At high-z, when the dominant 
star-forming galaxy is the central object and halos generally 
only contain a single large galaxy (i.e. the h alo occupation 
distribution lies on the plateau portion; e.g. iBerlind et al. I 
this is a reasonable assumption. At low-2, as halos 
grow, satellites become more numerous and house more of 
the total star formation, this becomes a poorer assumption. 
A galax;y driving an outflow does so locally, and does not 
necessarily know about the halo in which it resides. Hence 
in OD08 we moved to identifying individual galaxies and 
using their properties rather than the global halo proper- 
ties to derive wind parameters. In OD08 we presented a 
detailed comparison of the resulting global and IGM enrich- 
ment properties between our implementations, and found 
that at high-2 they gave similar results, but at low-z they 
differed. We believe our new implementation is more accu- 
rate particularly at z = 0, so we employ it here, and do not 
compare to our previous (OD06) implementation. 

A key aspect for studies of the ICM is iron en- 
richment from Type la supernovae, included in our 
new simulations. We employ a Type la rate following 
IScannapieco fc Bildsten" (l2005h. who parameterize observa- 
tions by iMannucci et ahl ( 2005t ) via a prompt component 
that tracks star formation plus a delayed component that 
tracks stellar mass. The formulae and metal yields are de- 
tailed in OD08. Metals and energy from delayed Type la 
SNe (coming from stars) are added to the nearest three gas 
particles. As noted in OD08, the vast majority of iron in the 
universe is produced in Type II supernovae, but for concen- 
trated stellar systems like groups and clusters. Type la's 
provide a larg e contribution to the gas-phase iron abun- 
dance (se e also Tornatore et al."1l2007l). 

As in lSpringel fc Hernguist] l|2003al ). each particle that 
is eligible for star formation h as its star formation rate com- 
puted based on a Kennicutt 

] 11998) law, and an appropri- 
ate level of enrichment is added based on Type II supernova 
yields (see OD08). The star-forming particle has some prob- 
ability to be kicked into an outflow or spawn a star particle. 
The ratio of those probabilities is given by the mass loading 
factor 77, which is an input parameter to the model. If a parti- 
cle is selected to be in an outflow, it is kicked with a velocity 
in a direction given by ±vxa (velocity cross accelera- 
tion), resulting in a quasi-bipolar outflo-s\[j. Wind particles 
propagate purely gravitationally until they are outside the 
star- forming region (typically few kpc), at which point they 
can again interact hydrodynamically with surrounding gas. 
The following formulae are used to relate the galaxy bary- 
onic mass Mb to the mass loading factor and wind speed for 



^ Movies can be seen at |http: //luca. as ■ arlzona.edu/~oppen/IGM/ 1 



star-forming particles in that galaxy: 



a = 200 



5 X IO12M0 n, 



Ho 



1/3 



kms 



4.3a V/l - 1, 
150km/s 



(1) 
(2) 
(3) 



where the first equation is taken from 
Murray. Quatert. fc Thompson I l|2005l ). originally based on 



Mo. Mao, fc White I (| 19981 ). T he value of /l is t aken from 
observ ations of outflows by iRupke. Veilleux"fc Sandore] 
l|2005h . who found that they lie between 1.05 — 2; we choose a 
random number in that interval for every outflow event. We 
include a factor to account for increased photon production 
in lower-metallicity systems. This wind mod el is identical in 
form to the "vzw" model of lOppenheimer' fc Dave (2006'J); 
full details are available in OD08. 

Ejected wind particles are hydrodynamically decoupled 
(i.e. the hydro forces are turned off, but not gravity) until 
they reach a density one-tenth of the critical density for star 
formation, or for a maximum period of 20/i^^kpc/i;u,, which 
is e.g. ~ 50 Myr for = 500 km/s. This is intended to 
simulate the outflow blowing a hole through the interstellar 
medium, which cannot be accurately represented given the 
smoothed nature of SPH. iDalla Vecchia fc Schavel (|2008| ) 
have investigated the effects of turning on or off hydrody- 
namic decoupling, and found that it gives significantly dif- 
ferent results for high-resolution simulations of disk galax- 
ies. However, we have done similar tests at the cosmologi- 
cal resolution used here, and found little difference in the 
global enrichment properties; this is perhaps not surpris- 
ing since the spatial resolution here is 5/i~^kpc, and the 
maximum distance an outflow particle can propagate dur- 
ing the decoupled phase is ~ 20/i~^kpc (and typically it is 
much less). One advantage of d ecoupling, as pointed out by 
iDalla Vecchia fc Schavel (|2008h . is that it results in better 
resolution convergence. Of course this may not be desireable 
if the result converges to a physically incorrect answer, but 
it at least allows a resolution-independent calibration of the 
outflow model. 

We consider two main runs in this paper. They 
are identical in every way (including random phases in 
the initial conditions) with the exception that one in- 
cludes momentum-driven outflows as described above, 
while the other does not include outflows. Note that in 
both cases, thermal energy feedback from supernovae is 
still accounted for in the m ulti-phase subgrid model of 
ISpringel fc Hernguist I l|2003ah ; however, the explicitly cou- 
pled phases means that such feedback cannot drive an out- 
flow. Thi s is the reason why an expli cit outflow model is 
required (jSpringel fc Hernguist Il2003bl ). 

Each run heis 256^ dark matter and 256'' gas particles in 
a randomly-generated volume o f (64ft~^Mpc)^ (como ving), 
with a concordance cosmology (|Komatsu et al. ||2008| ) hav- 
ing n = 0.25, A = 0.75, n = 0.95, Ho = 70 kms"^ Mpc'S 
erg = 0.83, and fit = 0.044. The value of ag is slightly higher 
than that favored from WMAP-5, moving more towards 
values derived from cluster abundances (e.g. lEvrard et al. I 
l2008l l: the difference is not important for the present work. 
Initial conditions are generated using an lEisenstein fc Hu] 
l|l999h power spectrum in the linear regime at z = 129 by 
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displacing particles off a grid according to the Zel'dovich ap- 
proximation. Particle masses are 2.72 X 1O*M0 and 1.27 x 
10^ Mq for gas and dark matter, respectively, and the grav- 
itational softening length is set to 5/i~^kpc comoving. 



2.2 Computing group properties 

To identify galaxy groups, we find large bound halos us- 
ing a spherical ove rdensity (SO) algorithm, as described in 
iKeres et al. I l|2005h . This algorithm finds particles at local 
potential minima, then expands spherically around those 
particles until the mean overdensity enclosed corresponds 
to the virial overdensity for t he assumed cosmology (specifi - 
cally, 178f7° *'5pcrit = 950nrit: lEke. Navarro, fc Frenk Ill998h . 
The outskirts of these systems may not properly follow the 
ellipticity of real systems (which is often ambiguous any- 
way) , but since the outermost properties contribute little to 
global group properties particularly when weighted by X-ray 
emission, this difference is not likely to be important. 

To identify gala xies, we use Splin e Kernel Interpola- 
tive Denmax fSKID: lKeres et al. Ilioosi '). Groups of galaxies 
are taken to be SO halos that contain ^ 5 resolved galax- 
ies. Our galaxy mass resolution limit is given by a stellar 
mass of M, > 8.72 x 10^ (i.e. 64 star particle masses), 
above whic h galaxy stellar prope rties have been shown to be 
converged (iFinlator et al. 2006) . For comp arison, this mass 
is « 0.2M„*„. (|Baldrv. G lazcbr ook. fc Driver 2008 ), where 



Afstar is the characteristic mass in a Schechter function fit 
to the present-day stellar mass function. 

The X-ray luminosity of individual gas particles is 
computed using t he Astrophysical Pl asma Emission CodsQ 
(APEC) models (|Smith et al. I I2OOII ). from the particle's 
mass, SPH-weighted density, temperature, and metallicity. 
APEC outputs X-ray spectra, from which the luminosity 
is computed by summing intensities from 0.5-10 keV, uni- 
formly weighted (i.e. no filter profile is assumed). Restrict- 
ing to 0.5-2 keV makes little difference for our systems. 
Contributions to line and continuum emission from individ- 
ually tracked elements (iron, oxygen, silicon, and carbon, 
along with hydrogen and helium) are computed separately 
and summed. Other elements are tied to iron a ssuming so- 
lar abundance ratios (| Anders fc Grevesse II 19891 ). though be- 
cause line cooling is dominated by iron and oxygen in the 
intragroup gas temperature and density regime, this choice 
makes little difference. 

The X-ray luminosity Lx of a given group is the 
summed X-ray luminosities of all particles within the virial 
radius. The group temperature Tx is computed as the X- 
ray luminosity-weighted temperature of gas particles, and all 
quoted metallicities are also X-ray emission- weighted. The 
emission weighting is done to mimic observations, which ob- 
tain temperatures and abundances by fitting APEC or sim- 
ilar models to X-ray spectra, and hence are effectively mea- 
suring luminosity-weighted quantities. Radial profiles are 
computed by performing similar summations in 25 equal 
thickness spherical shells about the SO halo center. 
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2.3 Resolution convergence 

In this paper we will not directly present any numerical res- 
olution convergence tests. This is unfortunate, and is not 
our preferred mode of operation. However, several practical 
issues limit our ability to directly study convergence for the 
objects of interest, primarily because the groups we exam- 
ine in this work are already at the extreme of what we can 
compute in a random-volume simulations. While we have 
smaller-volume simulations available, they do not produce 
any systems in the observed group-mass regime, so direct 
convergence tests on group properties are not possible. To 
be specific, our 64/i-iMpc run's largest group is 2 keV, 
which is already fairly small by X-ray group standards; our 
32/i-^Mpc run's largest is around 0.5 keV, which then does 
not allow us to check resolution effects for the observation- 
ally relevant 0.5 — 2 keV regime. On the flip side, we cannot 
run even larger volume simulations given our present compu- 
tational capabilities, because worsening the mass resolution 
beyond our 64ft^^Mpc run is found to significantly com- 
promise early galaxy formation, which is important because 
these large bound systems form much of their galaxy pop- 
ulation early on. Hence, we are stuck in the uncomfortable 
position of only having a single combination of mass reso- 
lution and volume that can examine the properties relevant 
for observed galaxy groups. For these reasons, the results in 
this work must be considered preliminary, and our main goal 
is to present new physical interpretations for group scaling 
relations as plausible if not definitive. 

Nevertheless, we do have some reason to believe that 
our results are robust. In OD08 we present resolution con- 
vergence studies for cosmic mass, energy, and metal distribu- 
tion, and the simulation resolution used here is reasonably 
well converged. We have also checked using a 32ft"^Mpc 
simulation with 256^ dark matter and gas particles that the 
group properties at temperatures below 0.5 keV are numer- 
ically converged. We hope to be able to run constrained 
simulations to form individual groups and clusters at higher 
resolution in the future, but for now we do not take that 
approach because we would like to compare large statisti- 
cal samples to observations. However, su ch a constrained- 
realization study for clusters was done bv iTornatore et al. I 
(2007) also using Gadget-2, with mass resolution compa- 
rable to what we use here, and through direct tests they 
determined that their results are not very sensitive to nu- 
merical resolution, at least in the central regions of clusters 
that dominate the X-ray emission. Hence we believe that 
our results here should be sufficiently robust to numerical 
effects, so that our conclusions should be generally correct. 



3 METAL CONTENT OF INTRAGROUP GAS 

Observations of cluster and group enrichment have been 
bolstered by the increased sensitivity provided by Chan- 
dra and XMM. It is well known that clusters show an 
iron abundance of one-third to one-half solar (outside the 
cool cor e ), ess e ntially independ e nt of cluster size (e.g. 
Renzini I Il997l : IPeterson etal"! l2003l : Ide Grandi et al. I 
2004h . Measurements in the poor group regime {Tx ^ 
2 keV), which have recently become available, show 
a larger scatter with a tendency towards slightly in- 
creased metallicity relative to rich clusters, although the 
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Figure 1. Top: [Fe/H]© (Lx-weighted) for galaxy groups drawn 
from our momentum-driven simulation (green) and no-wind sim- 
ulation (blue) at 2 = 0. Small numbers indicate the num- 
ber of resolved galaxies in the group. Data points (red) from 
iHelsdon fc Ponman I bOOd!) are shown for group s. The observed 
clust er iron abundance range of 0.33 — 0.7 solar l lPeterson et al. ] 
is indicated as horizontal dotted lines down to Tx = 1 keV. 
Middle: [Fe/H]Q in the momentum-driven wind simulation at 
2 = (green, reproduced from above panel) and 2 = 1 (magenta). 
The observed range at 2 = 1 for comparison is schematically 
indicated by the dotted lines showing 0.1 — 0.25 solar, i.e. half 
of the 2 = abundances. Bottom: [0/Fo]q in the momentum- 
driven (green) and no-wind (blue) simulations at 2 = 0, and 
the momentum-driven run at 2 = 1 (magenta). Range of [O/Fe] 
obser vations for clusters a t 2 Ri is indicated by the dotted 
lines llPeterson et al. II2003I) . 



poorest groups se em to reverse the trend and have 
lower abundances (Davis. Mulchaev &l Mushotzkv ] I1999I : 



iHelsdon fc Ponman I 200Cf ). It is puzzling that despite a sig- 
nificant increase in the ste llar baryon fraction from the rich 
clusters to groups (e.g. iGo nzalez. Zaritskv. fc Zabludoffi 
I2OO7I ). no corresponding change is seen in the X-ray gas 
metallicity. This suggests that the intragroup/intracluster 
medium metal enrichment is not simply related to the stars 
formed within the halo potential, and that distribution 
mechanisms and sources of enrichment play a critical role 
in governing metallicities. Hence the enrichment level of in- 
tragroup gas provides key cons traints on such processes (see 
review bv lBorgani et al71l2008l ). In this section we compare 
our simulations with and without outflows to observations 
of intragroup gas metallicities. 



3.1 Iron abundance 

Figure[T] (top panel) shows the iron abundance, weighted by 
X-ray luminosity, as a function of Lx-weighted temperature 
Tx for simulated groups. The green points show the results 
for our momentum-driven wind simulation, and the blue 
points show the no- wind run. The number of resolved galax- 
ies in each group is indicated by the small number on the plot 
(in this and subsequent figures). The key result is that the 
iron abundance hovers around the observed value between 
one-third and two-thirds solar for Tx ~ 0.5 — 2 keV sys- 
tems. The range of observat ions for clusters from 1 .5-10 keV 
([Fe/H]0 = 0.33-0.7) from lPeterson et al. I l|2003l ) is shown 
as the dotte d lines. Individual iron abu ndances for galaxy 
groups from IHelsdon fc Ponm"an1 (|200G|) are shown as the 
red po ints; observations bv lPavis. Mulchaev fc Mushotzkv I 
l| 19991 ) show a similar trend. There is a mild trend for lower- 
mass systems to have higher metallicities down to 0.5 keV, 
then a significant drop to smaller systems. Qualitatively, this 
is in agreement with observations, though data on the poor- 
est systems is still highly uncertain. 

Matching the observed enrichment level of intragroup 
gas is a significant success for this outfiow model, and in- 
dicates that our chemical enrichment model is working to- 
gether with stellar mass grow t h to properly enrich X-ray 
emitting gas. iTornatore et al.1 (|2007l ) achieved similar suc- 
cess for clusters, with a similarly sophisticated chemical en- 
richment model implemented in Gadget-2. Hence a metal- 
licity of one-third to one-half solar in intracluster and intra- 
group gas is achievable using standard yields and initial mass 
functions, when the cosmic production and distribution of 
metals is properly incorporated into hierarchical structure 
formation models. This suggests th at non-standard meta l 
yields or initial mass functions (IMFs: |Portinari et al. |[2004l ') 
are not required to mat ch o bserved abunda i ice lev els. In fact, 
ITornatore et al. I l|2007tl and lFabian et al. I l|2008h found that 
standard yields and a Salpeter IMF provides a better fit to 
cluster abundances and its evolution than a highly top-heavy 
IMF. 

Perhaps more surprising is that no wind simulation also 
roughly reproduces the observed iron abundance. One might 
surmise that this is because the iron content in intragroup 
gas is dominated by delayed Type la supernovae, thereby 
tracking the stellar mass distribution, and gaseous outflows 
play a subdominant role. To some extent that is true, but 
it turns out that the apparent insensitivity of iron enrich- 
ment to outflows is something of a coincidence. In reality, 
the no- wind simulation only matches the iron abundance be- 
cause it overproduces the amount of stars in groups, thereby 
compensating for the lack of metal dissemination into intra- 
group gas. We will return to this point when we examine 
the baryon fractions in i]3.3l 

The evolution of the iron abundance has now been mea- 
sured for massive clusters o u t to z ~ 1. iBalestra et al. I 
(200'A) and iMaughan et al. I (|2008l 'l found an increase of 
roughly a factor of two from 2 « 1 — > 0. The middle panel 
of Figure [T] shows the iron abundance in our momentum- 
driven case at 2; = (gree n points) and z = l(magenta). 
The dotted lines show the IPeterson et al. I l|2003l ) range at 
2 = lowered by a factor of two, for a rough comparison to 
the 2 = 1 predictions. Our simulation predicts iron abun- 
dance growth that is roughly consistent with that seen in 
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clusters, though the growth rate seems somewhat higher in 
the models than in the data. Of course we are comparing 
simulated groups with observed rich clusters at 2; = 1, so 
the discrepancy may be due to differential metal evolution 
between groups than clusters. But the predicted metal in- 
crease is independent of system size, at least in the range 
we can probe with these simulations, so this explanation 
is disfavored. Another possibility is that our central group 
galaxies have too much star formation at z < 1, since we do 
not include any feedback mechanism for truncation of star 
formation. These minor discrepancies aside, the predicted 
metallicity evolution is generally in agreement with obser- 
vations. 

The outflow simulation predicts a significant drop in 
metalhcity below 0.5 keV or so. This is likely a reflec- 
tion of the m ass-metallicity relation for galaxies setting 
in (e.g. iFinlat or & Dave 2008). Poor groups of this size 
typically house ^ L* galaxies (such as our Local Group), 
and below L« the metallicity of a galax y begins to drop 
with stellar mass ([Tremonti et al. II2004I ). Since the mass- 
metallicity relation is governed by outflows, one would ex- 
pect the no-wind simulation to show no such drop, and 
Figure [1] conflrms this to be the case. Hence although the 
star formation efficiency and stellar fractions are higher 
in smaller groups (as we will see in i]3.3p . outflows are 
able to compensate for this by driving more metals out 
of gr oup potentials. iDavis. Mulchaev fc Mushotzkv! (Il999l ) 



and iHelsdon fc Ponman I l|2000l ') found that the poorest 
groups in their sample at Tx ^ 1 keV have lower metal- 
licities (see data points at T < 1 keV in Figure [TJ, but the 
samples are small and the uncertainties large. Our simula- 
tion also predicts an increase in the spread in metallicities 
at low-Tx, owing to the stochastic nature of outflows; this 
may help explain the large range of metallicities seen for the 
poorest groups. As observations improve in this regime, the 
abundances of the smallest X-ray groups should provide a 
critical test of outflow models. 

3.2 Oxygen abundance 

The oxygen to iron ratio is often employed as a diagnostic 
of the contribution from Type II versus Type la supernovae, 
since oxygen is produced primarily in the former while iron 
primarily in the latter (at least in groups and clusters). The 
bottom panel of Figure [T] shows [O/Fe] for the momentum- 
driven wind simulation a.t z = Q (green points) and 2 = 1 
(magenta). These metallicities are X-ray emission weighted, 
so they trace the enrichment level of hot gas only. Since 
oxygen enrichment occurs predominantly at earlier epochs 
coincident with star formation, while Type la supernovae 
enrich in iron over longer cosmic timescales, [O/Fe] drops 
mildly with redshift. Our wind simulations predicts a ~ 30% 
drop from 2 = 1^0, independent of system size. Future 
observations of ICM alpha element abundances in high-z; 
clusters can test this prediction. 

Comparing [O/Fe] between our wind (green points) and 
no- wind (blue) simulations, it is evident that the [O/Fe] ra- 
tio provides a strong constraint on outflow properties. With- 
out outflows, oxygen is not disseminated effectively into dif- 
fuse gas, and instead remain trapped in stars. In our no- wind 
run, this results in [O/Fe]© ~ — 1 at all gr oup masses. Ob- 
servations of clusters from 1.5 — 10 keV bv lPeterson et alTI 



(2003!) indicate an oxygen to iron ratio of 70 ± 20% com- 
pared to solar, with no trend versus system size; this range 
is indicated by the dotted lines in the bottom panel of 
Figure [T] Our momentum-driven wind run places much 
more oxygen into X-ray emitting gas, and no trend with 
system size, in remarkably good agreement with observa- 
tions. It is notab l e that a different wind model explored by 
iTornatore et al." (12007 ), namely the "constant wind" model 
of lSpringel fc Hernquistl (|2003bl ). was not able to achieve as 
high [O/Fe] ratios as observed. Hence, [O/Fe] represents a 
sensitive test of outflows, which our momentum-driven wind 
simulation nicely passes. 

Contrary to the usual assumption, we flnd that the 
[O /Fe] ratio in X-ray emitting gas is not a good indicator of 
the relative numbers of Type II and Type la supernovae that 
have exploded within the group potential. Rather, [O/Fe] 
(and presumably other alpha-to-iron ratios) reflect more the 
distributions mechanisms of the various elements. This is 
seen by comparing the bottom panels of Figure [1] and Fig- 
ure [2I The latter shows the global mass- weighted [O/Fe] 
ratio in our simulated groups, including both stars and cold 
gas. In both the no-wind and momentum-driven wind sim- 
ulations, the mass- weighted [O/Fe] is around solar, roughly 
independent of group size. This is in stark contrast to typi- 
cal values of 0.7 solar in the outflow run, and 0.1 solar in the 
no-wind case. This arises because oxygen is predominantly 
distributed into hot gas via winds from galaxies, while iron 
is predominantly distributed from old stars. When no winds 
are present, oxygen remains trapped within galaxies, even- 
tually being locked into stars, while iron comes from older 
stars that can dynamically diffuse out of galaxies directly 
into hot gas. As a result. X-ray weighted [O/Fe] is a sen- 
sitive tracer of metal distribution mechanisms, but is not 
effective as a stellar chronometer. 

At face value, our results further imply that tidal and 
ram pressure stripping of metals out of galaxies falling into 
deep potential wells is not sufficient to pollute intragroup gas 
as observed. Such dynamical processes are in principle im- 
plicitly included in the no- wind run, yet the no- wind run fails 
to enrich intragroup gas sufficiently with oxygen. This result 
apparently co ntradicts the semi-an alytic-|-N-body simula- 
tion results of lKapferer et al~l (|2007l ), whose models predict 
that ram pressure stripping dominates hot gas enrichment 
for present-day clusters. This may indicate that the simpli- 
fied modeling used in that work is not capturing the full 
hydrodynamics. Alternatively, it may be that ram-pressure 
stripping is numerically supp ressed in our simulations. For 
instance, lAgertz et al. l|2007t ) demonstrated that Gadget- 
2 tends to oversuppress stripping relative to AMR codes, 
as SPH has trouble modeling hydrodynamical instabilities. 
Hence while our model of dynamical metal ejection is suc- 
cessful, further testing is required to determine whether gas 
stripping processes can also be effective. 

3.3 Baryon fractions 

Figure [5] shows the baryon fractions as a function of halo 
mass out to rsoo in our simulations with winds (green points) 
and without (blue). The top panel shows the fraction of 
baryons in stars (we do not distinguish stars in galaxies from 
intracluster/intragroup stars). It is evident that galactic out- 
flows are effective at suppressing star formation on group 
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Figure 2. Baryon fractions is various phases, out to rsoo- Green 
points are for our favored momentum-driven wind simulation, 
while blue points are for the no-wind run. Top: Stellar baryon 
fraction, including all stars in the groups regardless of whether 
they are within galaxies or not. Red points and red line show 
data and fit to obs ervations of stellar baryon fractions from 
iGonzalez. Zaritskv. fc Zabludoff (2007); their data for groups 
< 10^'* M0 are not shown because they lie beyond the range where 
the M500 — rsoo relation has been calibrated, so are highly un- 
certain (A. Gonzalez, priv. comm.). Top axis shows approximate 
X-ray temperature based on a fit to A/500 ~ Tjf from our sim- 
ulations. Second: Cold baryon fraction, where cold gas is all gas 
below lO^K. Third: Baryon fraction with respect to total mass. 
Dotted line shows the ratio Hj/Om = 0.176 assumed in our sim- 
ulation. Bottom: Log of oxygen to iron mass to rsoo, including 
both gas and stars, relative to solar ratio. This can be compared 
to the Lx -weighted [O/Fo] in Figure [Tl 



contributions (|Baldrv. Glazebrook. fc Driver I [iooi ) . Below 
this mass, systematic uncertainties become large. 

The most massive groups still have 20% of their 
baryons in stars, which is too high compared to data, 
and a simple extrapolation of the simulated group trends 
to higher masses would clearly violate observations of 
stellar fractions in rich clusters (typically <g 10%). An- 
other way to say this is that the observed /* — Afsoo re- 
lation (where /* is the stellar fraction) is steeper than 
our simulations predict. Specifically, our simulations yield 
d log /* /d log M500 ^ 0.18 around 10^"^ Mp, which is shal - 
lowe r than both [Gonzalez. Zaritskv. fc Zabludoff I (|2007l ) 
and iLin fc Mohr I ( 20041 ) data, though some what steeper 



than the semi-analytic mod el predictions of iBower et al. 



.2006) (as calculated by iBaldrv. Glazebrook. fc Driver 



20081 . However, this does not necessarily indicate a sig- 
ni ficant constraint on cold dark matte r models (as argued 
bv lBaldrv. Glazebrook. fc Driver |2008| ) . Rather, it probably 
owes to a lack of suppression of star formation in massive 
systems in simulations, from AGN feedback or whatever the 
heating source is that prevents cooling flows. 

The suppression of star formation on group scales is 
not as great as the global suppression: The simulation with- 
out winds at z = contains 17% of its total baryons in 
stars, while the wind run has 7%. The latter value is in 
good agre ement with observational est imates (|Bell et al. I 
I2OO3 ; Baldrv. Glazebrook. fc Driver I liooa ). Hence outflows 
are able to suppress star formation both globally and on 
group scales in general agreement with observations. 

The third panel of Figure [2] shows the mass frac- 
tion of total baryons against halo mass. In general, the 
total baryon fraction is slightly smaller than the cos- 
mic baryon fraction (indicated as the dotted line), ow- 
ing to increased pressure support that puffs out the hot 
gas; this is qualitatively consistent wit h observations (e.g. 
iGonzalez. Zaritskv. fc Zabludoff I |20o7). Outflows are fur- 
ther able to remove a small but noticeable amount of baryons 
from halos with masses < IO'^^Mq. 

Returning to the iron abundance with and without out- 
flows, it is now clear that although simulations without out- 
flows reproduce the observed intragroup metallicities, they 
do so by forming too many stellar baryons. Outflows there- 
fore serve the dual purpose of suppressing star formation, 
while transferring sufficient metals from galaxies into the 
hot gas phase. Furthermore, outflows significantly increase 
the alpha element ratio in intragroup gas to be in better 
agreement with observations. Hence our simulations suggest 
that outflows are likely to be necessary in order to enrich 
intragroup gas to observed levels. 



mass scales. Observational estimates of the stellar fraction 
on group scales are challenging, primarily owing to the dif- 
ficulty in constraining th e dominant baryonic component 
in hot gas, but data from IGonzalez. Zaritskv. fc Zabludoff 1 
l|2007l ). which carefully accounts for intracluster stars, are 
broadly reproduced in the momementum-driven wind run. 
As seen in the second panel, there is virtually no cold gas 
(T < lO^K) at least in the range of observed X-ray groups, so 
baryons not in stars are essentially all in an X-ray e mitting 
plasma. The lGonzalez. Zaritskv. fc Zabludoff 1 1 20071) stellar 
fractio n data are consistent with data from iLin fc Mohr I 
l|2004l ) for M500 > 10" Mq when corrected for ICL light 



4 PRE-HEATING 

We now turn to studying another aspect of outflows, namely 
the energy they add to intragroup gas. Entropy levels in 
intragroup gas are known to be elevated relative to self- 
similar extrapolations from clusters, but th e physical pro- 
cess r esponsible for pre-heating is unclear. IVoit fc Brvan I 
(|200ll ) argued analytically that supernova heating together 
with radiative cooling could yield the observed entropy dis- 
tribution in groups. Incorporating such feedback processes 
into simulations has only recently been attempted; for in- 
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Figure 3. Entropy, defined as Tx/ne^^i as a function of X-ray 
temperature for galaxy groups in our simulation with outflows 
(green) and without (blue). Binned data from PSF03 is shown in 
red. Expected self-similar scaling, normalized to the most massive 
data points, is shown as the magenta line. Top: Entropy at 10% 
of the virial radius. Bottom: Entropy at rsoo, roughly two-thirds 
the virial radius. 

stance, iBorgani et aTl feOOS*) examined and dismissed a sim- 
ple mode l of galactic outflows as the dominant source of pre- 
heating. iMuanwong. Kay, fc Thomas I l|2006l ) showed that 
an ad hoc model of energy input proportional to star forma- 
tion is able to reproduce observed entropy levels, but pre- 
liminary comparisons to entropy evolution favored an early 
pre-heating scenario. None of these works employed an out- 
flow model that is constrained to match external observa- 
tions. Here we examine whether our momentum-driven wind 
model, in addition to enriching intragroup gas, is concur- 
rently able to add significantly to its entropy budget. 

4.1 Entropy 

Figure [3] (top panel) shows entropy 5o.i at O.liZvir as a 
function of group X-ray temperature, for our momentum- 
driven wind (green) and no-wind (blue) runs. Entropy 



is defined here as Tx/n^^^, where Tx is the system X- 
ray luminosity-weighted temperature and rie is the elec- 
tron density at a given radius. Data points (binned) from 
PSF03 are shown, and the expected self-similar scaling nor- 
malized to the most massive clusters is indicated by the 
magenta line. In self-similar scaling, rie (O.liivir) is inde- 
pendent of size since groups are scaled-down versions of 
clusters, hence S cx Tx- As is well-known, So.i is ele- 
vated at group scales relative to a self-similar extrapo- 
lation from ma ssive systems. While an entropy fl oor at 
~ 100 keV cm^ (|Ponman. Cannon, fc Navarro II 19991 ) is not 
seen in the models, more recent observations indicate no 
floor but merely a more gradual decline than expected. 
PSF03 finds a best fit of approximately So.i oc . 

Elevated entropies are obtained in both wind and no- 
wind runs, since radiative cooling is effective at preferentially 
rem oving Iqw-entropy intragroup gas over a Hubble time 
(e.g. lBryanl2000l : IVoit fc Bryan I2OO1I) . The low-entropy gas 
then ends up forming stars. While the entropy trend pro- 
duced by r adiative cooling qualitatively agr ees with data, 
as shown in lDave. Katz. fc Weinberg ' ("2002!), there are two 
major difficulties with the radiative cooling solution alone. 
The first is that the amount of stars produced by the low- 
entropy gas exceeds o bservations, as shown in Figure [2] and 
found previously (e.g. iMuanwong et al~ll2002l ). The second 
is that, as seen in Figure |3l the entropy increase from radia- 
tive cooling alone is not quite sufficient to match observa- 
tions. Hence an alternate source of entropy is still required, 
preferably one that simultaneously reduces the amount of 
star formation in groups. 

Our momentum-driven outflow model appears to allevi- 
ate both problems. We have seen that it reduces the stellar 
baryon fraction to be in better agreement with observations 
on group scales. Figure |3] also shows that it raises the en- 
tropy level over the no- wind case, particularly for systems 
with Tx < 1 keV, to be in slighly better agreement with 
the low-Tx end of available observations. An extrapolation 
of the observed trend (5* oc T"^^) to lower Tx is also well- 
matched to the outflow run predictions, which can be fit 
by 5*0.1 cx: Tx^ over the temperature range probed. Hence 
outflows together with radiative cooling appear to provide 
a viable solution to the entropy crisis in galaxy groups. 

The fact that outflows increase core entropy is not triv- 
ial, since they also suppress cooling, and hence suppress the 
removal of low-entropy gas. The entropy addition from out- 
flows therefore appears to be (more than) compensating for 
the reduced stellar fraction. That said, our most massive 
systems show twice the stellar fraction as observed (see Fig- 
ure 121), and so some additional mechanism such as AGN 
feedback is likely required to further suppress star forma- 
tion at Tx ^ 1 keV. It is not clear whether such a mecha- 
nism would continue to compensate for further reductions in 
stellar fractions in the way that outflows do. Looking care- 
fully at Figure 131 it does appear that the core entropy is still 
slightly low for Tx > 1 keV systems, though it is within 
uncertainties. Hence it could be that radiative cooling, out- 
flows, and some other feedback mechanism all work together 
to establish the entropy level in larger groups and clusters. A 
fully self-consistent model incorporating all these processes 
is required to assess this. 

The entropy in group outskirts is also seen to be ele- 
vated with respect to a self-similar extrapolation from clus- 
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ters. The bottom panel of Figure |3] shows the entropy at r^oo, 
roughly two-thirds the virial radius. Observations (PSF03, 
shown as the red points) are more uncertain here, partic- 
ularly for groups, but our simulations generally fall in the 
observed range. Saoo is mostly insensitive to outflows; it is 
apparently established by radiative cooling. If anything, the 
wind run shows marginally lower entropies, probably owing 
to reduced removal of low-entropy baryons. The S500 — Tx 
relation shows a slope just below unity, unlike 5*0.1 and more 
in line with expectations from self-similarity. If this would 
continue to the largest clusters, it would exceed the observed 
Ssoo; this will require constrained realization simulations of 
large clusters to test. 

4.2 Outflovir energy deposition 

We now quantify pre-heating due to outflow energy injec- 
tion in simulations. During our runs we record the mass and 
velocity (and hence energy) of particles kicked into an out- 
flow. Ideally, one would like to track where and when this 
energy gets deposited into surrounding gas; unfortunately, 
this is not straightforward. So instead, we make a crude es- 
timate of energy input into intragroup gas by assuming that 
all the outflow particles lying within a group's hot gas by 
2 = have deposited their entire outflow energy into that 
gas. Note that this estimate may either be an underestimate 
of the actual energy added if outflow particles lose energy to 
intragroup gas but end up outside the group, or it may be an 
overestimate if groups accrete outflow particles that already 
have deposited most of their energy outside the group's halo. 
While this estimate is not exact, it is still instructive. 

Figure |4] shows the wind energy per hot baryon de- 
posited as a function of Tx (top panel), and the median 
age when the energy was deposited (bottom panel), as the 
green points. The energy is the total outflow energy of all 
wind particles residing within the hot intragroup gas, di- 
vided by the number of baryons in that gas (note that this 
is not all the baryons in the group). The age is estimated 
by taking the median of the ejection times of all wind par- 
ticles within intragroup gas; again, this must be considered 
a relatively crude estimate of the energy injection time. 

The energy deposition per baryon shows a mild decreas- 
ing trend with system size. This is at first glance surprising, 
since in the momentum-driven outflow model the wind en- 
ergy output per unit star formation scales as rjv^ oc Ucirc 
(equation|3]). Hence one expects larger systems to have larger 
energy input. But the baryonic budget serves to counter this 
trend: Larger groups have a lower stellar fraction and more 
hot gas over which to spread the energy (Figure [2]). The 
net result is that hot baryons in larger systems are heated 
slightly less by winds. 

The total energy injection into groups from outflows is 
comparable to the total supernova energy from stars formed 
in the group (magenta points in Figure 3}. The SN energy is 
computed by taking the total stellar mass in each group and 
multiplying by 4x 10** ergs/M© for a standard IMF (OD08). 
As argued in OD08, such high levels of energy input appear 
to be necessary to provide the cosmic distribution of met- 
als from galactic outflows, and owing to expected radiative 
losses this level of energy may not be obtain able from su - 
pernovae alone (barring IMF variations; e.g. iDave I [20081 ). 
While we cannot claim to fully understand or model the 
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Figure 4. Top: Estimated energy per baryon deposited into in- 
tragroup gas from galactic outflows (green points), as a function 
of system temperature. Solid line is the line of equality between 
wind energy and system thermal energy. Magenta points show 
supernova energy per hot baryon from all the stars formed within 
the groups. Bottom: Age by which half the energy was deposited 
(green points), taken to be the median ejection time of all wind 
particles in the hot intragroup gas. Magenta points show median 
age of stars in group. 



outflow driving mechanism, we note that such energetic out- 
flows seem to be required to match a wide range of observa- 
tions, not just intragroup gas entropy. 

The median energy injection time (green points in bot- 
tom panel of Figure |4)| is typically fairly recent, with the 
smallest systems having a median pre-heating age of 5 Gyr 
[z « 0.5), while our largest systems at around 2 keV have 
a typical pre-heating age of 8 Gyr [z ~ 1). This can be 
compared to the median age of star particles in the group 
(magenta points), which are comparable although slightly 
older, particularly in the smaller systems. In general it ap- 
pears that stars in group galaxies are responsible for en- 
riching and pre-heating intragroup gas, as expected. In the 
future we plan to examine the redshift evolution of hot gas 
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to determine more precisely where and when outflow energy 
and metals are deposited. 

4.3 Scaling relations 

Figure |4] shows that the energy deposition from winds typi- 
cally exceeds the system's thermal energy for systems below 
about 0.5 keV (note the solid line showing equality between 
Tx and wind energy input). Hence outflows should produce 
a lowered luminosity for systems near this temperature and 
below. Here we investigate whether this provides an expla- 
nation for the apparent break in luminosity scaling relations 
near ~ 1 keV, below which the Lx — Tx relation is observed 
to be significantly steeper (e.g. IXue fc Wu II2OO0I ). 

Figure [5] shows the Lx — Tx (top) and a — Tx (bottom) 
relations for our simulated groups from our momentum- 
driven wind run. The Lx — Tx relation shows a clear 
break around 1 keV, with the power-law slope above and 
below this break being 2.6 ± 0.6 and 4.9 ± 0.2, respec- 
tively. These slopes agree well with observations: The red 
dashed li ne shows a fi t to o bservations of groups and clus- 
te rs fromlXue fc Wu I ||2000|) ; a fit to Lx - Tx for clusters 
bv lWhite. Jones, fc Forman I l|l997l . not shown) is quite sim- 
ilar. Green points sh o w da ta from a ROSAT compilation by 
lOsmond fc PonmanI (|2004 l: we show all their groups that 
have measured Tx values, including systems with a small 
radial extent in observable X-rays (i.e. their "H-sample") 
that they exclude from their fits, since we make no such cuts 
in our simulations. The Lx values used are extrapolated to 
fsoo, whereas our simulations include all luminosity out to 
the virial radius, but this sho u ld not be a large difference. 

The lOsmond fc PonmanI ||2004 ) Lx - T - X data as 
well as the lXue fc Wu I l|2000l ') fit at Tx < 1 keV lie sys- 
tematically below that for the simulated groups. If real, this 
suggests that there may still be some missing energetic input 
into low-mass systems. 

The no-wind simulation (not shown) also displays a 
break in the Lx — Tx relation, but it is not as promi- 
nent and occurs at a lower Tx- Our no- wind case is in gen- 
eral agreement with simu l ations (also without outflows) by 
iDave. Katz. fc Weinberg] (|2002l ). who argued that cooling 
alone can produce a break in scaling relations, but this break 
occurs at ^ 0.3 keV. Hence outflows serve to add enough en- 
tropy to move the Lx break closer to the observed location. 

The cr—Tx relation (bottom panel of Figure[S} shows no 
break. Here, a is calculat ed from the dark matter particles in 
each simulated group; in lDave. Katz. fc Weinberg I (|2002| ) it 
was shown that calculating a from galaxies yields a similar 
values so long as there are a sufficient number of galaxies in 
a group. The simulated groups are best fit by a power-law 
slope of 0.60 ± 0.02 (solid line), and the fit is remarkably 
similar to that based on compiled observations of groups 
and clusters bv IXue fc Wu I (|200G| . red dashed line). This 
is slightly different from the self-similar scaling slope of 0.5, 
and shows that groups at a given mass have elevated temper- 
atures. Th e no-wind simulation (not shown ) shows a similar 
slope (see iDave. Katz. fc Weinberg I l2002h . indicating that 
the elevated temperatures occur primarily because radiative 
cooling has preferentially removed the low-temperature gas. 

The scatter in the observational points is considerably 
larger than in the simulated case. As an example, green 
points (with errors) on the Lx — Tx plot show individual 
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measurements from lOsmond fc PonmanI (|2004l ) . The large 
scatter may be due partly to observational systematic un- 
certainties, but this seems unlikely to provide the bulk of 
the explanation, since clusters show similar scatter that can 
be directly linked to their well-measured den sity and tem- 
perature profiles (e.g. iMarkevitch et al. |[l998l ). It has been 
suggested that it m ay instead reflect int rinsic scatter in the 
level of pre-heating ( Balogh et al. |[2006l ). Since our simula- 
tions implicitly produce variations in halo concentration and 
formation epoch, as well as in wind energy deposition, these 
alo ne cannot provide t he observed scatter (as also argued 
bv iBalogh et al.1 l2006il ). We leave for the future a careful 
consideration of observational causes of scatter via making 



12 R. Dave, B. D. Oppenheimer, S. Sivanandam 



mock X-ray images and analyzing them alongside real X-ray 
data. This will also be useful for evaluating the more general 
bias that observed groups tend to be X-ray bright, but these 
may not be representative of groups as a whole. The final 
solution, however, may require some intermittent additional 
source of heating such as AGN feedback. 

In summary, outflows add an important contribution 
to the entropy budget of poor groups, whose energy in- 
put becomes comparable to the system thermal energy at 
Tx ~ 1 keV. Radiative cooling is primarily responsible for 
breaking self-similarity above this scale, while outflows ex- 
aggerate deviations from self-similarity below this scale. 



5 PROFILES 

With increasingly sensitive observations from Chandra and 
XMM, it is now possible to construct radial profiles of group 
properties. Such profiles provide an added challenge to mod- 
els, one that simulations have traditionally had a difficult 
time matching, since they contain more detailed information 
about metal and energy injection processes. In this section 
we examine projected radial profiles for physical quantities 
from our momentum-driven wind simulation. 

Figure |6] shows the luminosity, temperature, electron 
density, entropy, iron abundance, and oxygen to iron ratio 
profiles as a function of radius, from « 0.06 — l-Rvir- The 
inner limit corresponds to ~ O.I-R500, which is typically few 
tens of kpc, and approaches the efi^ective resolution limit of 
our simulations. As seen in the top panel, the luminosity 
profile shows a flat ce ntral core (similar to observed non- 
cool core systems; see iDave. Katz. fc Weinberg Ibooj ). so 
the poorly resolved central region does not make a signiflcant 
contribution to groups' X-ray emission-weighted properties. 

We subdivide the groups into two equal-number bins 
exceeding 0.5 keV, which results in about 40 groups each 
from 0.5 - 0.7 keV and 0.7 - 2 keV. We then compute the 
median value among the 40 groups for each physical quantity 
within each of the 20 radial bins. This yields the median pro- 
flies shown within these two temperature ranges. We show 
such profiles for our momentum-driven wind run (green and 
magenta) and no- wind run (blue and cyan). 

The simulated temperature profiles show a drop by 
a factor of ~ 2.5 from 0.1 — > l-Rvir- The no-wind pro- 
files a re well-described by the fit proposed bv iLoken et al. I 
l|2002h . namely Tx oc (1 -f 0.75r/r5oo)"^ ^ With outfiows, 
the profiles are altered slightly in the inner parts, such that 
they show a flattening and perhaps even a turnover within 
~ 0.2_Rvir. Qualitativel y, this is in better agreement w ith X- 
ray group proflles from iRasmussen fc Ponmani (|2007l '). who 
flnd a turnover in the temperature proflle at ~ O.lfJsoo- 

Fits to observed temperature profiles (T/ < T >)are 
shown as the red and blue dotted lines in Figure |6l arbi- 
trarily scaled to lie in t he ran ge of our simulated groups. 
IRasmussen fc Ponman I (|2007l . red dashed line) measured 
profiles for 0.3 — 2.1 keV groups, so this is most di- 
rectly compa r able t o our simulated systems. We show the 
iLoken et al. I l|2002t) form of the fit to their data, specifi- 
cally Tx oc (1 -I- 0.75r/r5oo)^^ which is somewhat shal- 
lower than the si mulated profiles. We also s how a fit to 
cluster profiles bv iLeccardi fc Molendi I (|2008l . blue dashed 
line). Their analysis is notable because of the careful ac- 
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Figure 6. Radial profiles of physical quantities in our 
momentum-driven wind simulation. Profiles shown are medians 
over all galaxy groups within two bins: 0.7 < Tx < 2 keV 
(green) and 0.5 < Tx < 0.7 keV (magenta); each bin con- 
tains 40 groups. Blue and cyan curves show median profiles 
of groups in those two temperature ranges from the no- wind sim- 
ulation. Note that there is significant scatter in individual profile 
shapes, which is not shown. Dott ed lines in Tx and [Fe/H ] panel s 
shows observ ations; red are from IRasmussen &: Ponmani l|20o3), 
blue are from lLeccardi fc Molendi I l|2008l l. Toy: X -ray luminosity 
(erg/s); Second: X-ray luminosity- weighted temperature (keV); 
Third: Free electron density (cm~^); Fourth: Entropy (keV cm'^); 
Fifth: Luminosity-weighted iron abundance; Bottom: Oxygen to 
iron ratio. 



counting of various systematic effects such as background 
conta mination and biases from low surface brightness re- 
gions (iLeccardi fc Molendi II2007I). b ut other groups' results 
are similar (|de Grandi fc Molendi I I2OO2I : IVikhlinin et al. I 
2005). Although the observed clusters are not as directly 
comparable to our simulated groups, they may be more ac- 
curate particularly in the outer regions, and since the profile 
shapes seem to be mostly independent of system size, they 
may still be relevant. 

The simulated Tx profiles appear mildly steeper than 
observations at r > O.l Ruir, as has historically been seen in 
such comparisons, (e.g. iMarkevitch et al. |[l998l ). There may 
be some additional source of thermalization such as cosmic 
rays or conduction that serves to flatten the real temper- 
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ature profiles. So long as such a process only moves heat 
and not mass or metals around, then the remaining physi- 
cal profiles should not be significantly affected by fixing this 
discrepancy. At smaller radii, we do not have the resolution 
to reliably probe the temperature profile too far in, but the 
observed profiles show a drop towards the center that is sim- 
ilar to but steeper than seen in simulations (at least for the 
larger groups) . This may result from an excess of star forma- 
tion in the simulated central group galaxies are over-heating 
the core. Overall, the agreement with temperature profiles 
are reasonably good but not perfect. 

The observed iron abundance profiles are in reasonable 
agreement with observations. There is a hint that the sim- 
ulated abundance profiles are steeper than observed, but 
in this case the uncertainties are large enough (particularly 
in the outer regions) that the differences are probably not 
significant. Our simulations may also not be resolving the 
outer iron abundance profile (r > 0.37?vir) at our numerical 
resolution (jTornatore et aLlboOTl ). Outflows do not greatly 
impact the iron abundance profile, because the stellar distri- 
bution is not altered significantly by outflows even though 
they suppress star formation overall. 

The hot gas density profile shows a clear impact from 
outflows at Tx < 1 keV. The We profile is significantly shal- 
lower in the central region for lower-mass groups in the out- 
flow run. In contrast, the no- wind run shows little difference 
between the two mass bins, although the increased efficiency 
of low-entropy gas removal does cause the smaller systems 
to have slightly flatter inner profiles. This demonstrates that 
outflows are important for further breaking self-similarity in 
the core regions of poor groups as observed. 

The entropy profiles are also affected by outfiows. The 
profiles from the wind run are significantly shallower than 
those from the no-wind run. Furthermore, less massive 
groups show a shallower profile relative to more massive 
systems. The outer profiles generally approach a slope of 
1.1 expected for shock heating, but smaller systems do so 
further out. PSF03 measured entropy profiles that look qual- 
itatively similar, showing a slope shallower than 1.1 out to 
^500, particularly for smaller systems. 

Finally, the alpha enrichment profiles (here traced by 
[O/Fe]) show a clear difference in amplitude between no- 
wind and wind runs. This reflects the overall alpha enrich- 
ment difference seen in Figure [T] More subtly, the no- wind 
runs show almost no gradient in [O/Fe], while a mild gra- 
dient is seen in the wind run. Still, the gradient is generally 
not as strong as observed in groups. iRasmussen fc Ponman I 
|2003) found that silicon to iron increased by a factor of 
2.5 from 0.1 — l-Rsoo, whereas our predicted [O/Fe] only in- 
creases by 40% over a comparable range. This may be owing 
to excess central star formation in our simulated groups that 
raises the central oxygen abundance. 

Overall, the X-ray profiles give a more detailed view of 
the physics affected by outflow energy input as a function 
of mass. Winds are seen to clearly affect the gas density 
profile (and hence the entropy proflle), as well as the alpha 
element abundance proflle. Comparisons with observations 
yield general agreement, though in detail signiflcant differ- 
ences remain. Hence these simulations do not completely 
explain all observed group properties, even though global 
X-ray luminosity-weighted properties are fairly well repro- 
duced. As detailed X-ray profile observations of groups im- 



prove, matching them should yield insights into other feed- 
back or heat conduction processes operating within intra- 
group gas. 



6 CONCLUSIONS 

Using ACDM cosmological hydrodynamic simulations of 
galaxy formation including galactic outflows and a sophisti- 
cated chemical enrichment model, we have studied the distri- 
bution and evolution of metals and entropy in galaxy groups. 
The one-line summary is that our simulation with outflow 
scalings as expected for momentum-driven winds is broadly 
successful at reproducing both the enrichment and entropy 
level seen in X-ray emitting intragroup gas. In more detail, 
we find: 

• Our model with outfiows reproduces the observed iron 
abundance of intragroup gas, while maintaining a stellar 
baryon fraction in broad agreement with observed values. 
The iron abundance is somewhat higher than seen for rich 
clusters, but then drops for Tx < 0.5 keV systems. 

• Predictions for [O/Fe] in X-ray emitting gas are 
strongly affected by outflows. Our momentum-driven wind 
model is successful at reproducing the [O/Fe] ratio observed 
for clusters. Without outflows the oxygen abundance is too 
low, although it may be that our SPH-based code is un- 
derestimating the effectiveness of gas (and metal) stripping. 
Measurements of [O/Fe] in the poor group regime should 
provide stringent constraints on outflow models. 

• The mass-weighted and X-ray luminosity weighted 
[O/Fe] ratio within groups can be signiflcantly different. 
This suggests that [O/Fe] in X-ray luminous gas is a poor in- 
dicator of the relative contributions of Type II and Type la 
supernovae within the group as a whole, and is instead gov- 
erned primarily by the effectiveness of distribution mecha- 
nisms for oxygen and iron. 

• Iron abundance evolution predicted for groups is in gen- 
eral agreement with that seen for rich clusters, i.e. roughly 
a factor of two increase from z — 1 ^ 0. Our simula- 
tions predict that [O/Fe] decreases by ~ 30% over the same 
time span, though halting late-time star formation in central 
group galaxies may increase such evolution. 

• The total baryon fraction in groups is predicted to be 
somewhat below the cosmic mean baryon fraction, owing to 
pressure support from hot gas. In our outflow model this 
deficit increases in lower mass systems owing to ejection of 
halo baryons, while without winds the baryon fraction ex- 
ceeds the cosmic value in the poorest groups. 

• The core (O.l-Rvir) entropy level of poor groups is found 
to scale as 5o.i oc Ty^, in agreement with observations by 
iPonman. Sanderson, fc Finoguenov I l|2003l ) . Much of the en- 
tropy increase over self-similar scaling owes to radiative cool- 
ing alone, but outfiows add significant entropy to systems 
below 1 keV. 

• Entropy at -R500 is elevated from self-similar scalings, 
in broad agreement with observations. However, the scaling 
with temperature within the poor group regime is close to 
self-similarity. Outflows have essentially no impact here, as 
the elevation of S500 owes purely to the removal of low- 
entropy gas via radiative cooling. 

• The amount of energy per baryon added to hot intra- 
group gas drops with increasing system size, and exceeds 
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the thermal energy below 0.5 keV. As a result, the Lx — Tx 
relation shows a clear break below 1 keV. Both the Lx — Tx 
and a — Tx relations predicted from our momentum-driven 
wind run are in good agreement with observations. 

• Profiles of Tx and [Fe/H] are qualitatively similar to ob- 
served, though model profiles appear to be somewhat steeper 
than data. The alpha enhancement profile is also quite shal- 
low. Possible causes include missing heat transport mecha- 
nisms and an excess of late-time star formation in the central 
galaxy. 

• Outflow energy injection makes the density and entropy 
profiles shallower at Tx < 1 keV, particularly in the inner 
regions. 

The simultaneous success of enrichment and entropy in- 
jection from our momentum-driven wind outflows is a signif- 
icant achievement, especially considering that this outflow 
model was not tuned in any way to match observations of 
intragroup gas. While the iron abundance is not very sen- 
sitive to outflows, the oxygen to iron ratio provides a key 
diagnostic that is well reproduced by our outflow model. 
Most of the entropy increase actually owes to radiative cool- 
ing, but outflows serve both to enrich intragroup gas as well 
as lower the stellar content to bring it more into line with 
observations. 

The implication of our models is that pre-heating in in- 
tragroup gas arises from radiative cooling together with the 
same galactic outflows required to, for instance, enrich the 
IGM at high redshifts or establish the mass-metallicity re- 
lationship of galaxies. Hence the source of pre-heating may 
not be a mysterious widespread energy injection at high red- 
shifts, nor be related to AGN feedback at lower redshifts. 
While something like AGN feedback is still required to trun- 
cate star formation in the central galaxy, our results suggest 
that it need not be a significant contributor to the overall 
entropy budget of intragroup gas. This eases the tension be- 
tween a feedback process that must operate on dense central 
gas to prevent cooling flows (and therefore requires large 
amounts of energy to significantly increase entropy, since 
S oc n"^''^), and the need to raise the entropy of the more 
diffuse intragroup medium. 

Nevertheless, there are certainly still significant prob- 
lems with our simulations as-is, most notably the excess of 
late-time star formation in the central galaxy that may be 
responsible for poor fits to observed X-ray profiles. There is 
also the issue of the observed scatter in the X-ray scaling 
relations, which points to some more stochastic form of en- 
ergy injection. Finally, it appears that large galaxy clusters 
also h ave elevated entropies a ll the way out to their outer re- 
gions l|Dunn fc Fabian"ll2008l ). Since wind energy cannot be 
responsible because it is less important in higher mass sys- 
tems, this suggests that AGN energy, perhaps transported 
via bouyant bubbles, are needed to explain the thermody- 
namics of clusters. 

A complete understanding of hot diffuse gas in groups 
and clusters will only come when a model addresses all three 
IGM issues simultaneously, namely the cooling flow, entropy, 
and metallicity problems. Our current simulations explicitly 
do not do so, focusing on two of the three problems solely 
in the poor group regime. Hence although our simulations 
suggest that entropy injection is tied mostly to metal en- 
richment and not cooling flow prevention, we forward this 



interpretation only a suggestion and not a definitive state- 
ment. This must still be tested using simulations that self- 
consistently incorporate some mechanism to quench star for- 
mation in massive galaxies in accord with cluster galaxy 
data. We are actively pursuing such investigations. 

There are two immediate paths for future investiga- 
tions. The first is to incorporate some physical or heuristic 
model of star formation truncation in massive galaxies, in 
order to understand its impact on energy and metal input 
into intragroup gas. The second is to study how the en- 
richment and entropy of intragroup gas grow with redshift 
within a hierarchical framework, in order to better under- 
stand its physical origin. We also plan to continue developing 
tools to provide closer comparisons with observations, such 
as the ability to produce simulated X-ray maps that can be 
analyzed alongside data. Finally, we are planning to run sig- 
nificantly larger simulations with up to 2 x 5i2^ particles in 
the near future, which will allow us to have greater dynamic 
range, study larger clusters, and properly assess resolution 
effects. While our initial efforts towards understanding intra- 
group gas chemo- and thermodynamics from outflows seem 
promising, there is much work to be done to fully under- 
stand the interplay between galaxy formation processes and 
the X-ray emitting plasma in the largest bound structures 
in the Universe. 
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